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ABSTRACT 

The effect of solar rotation on the proton-alpha differential flow speed, Vap, and consequently on the 
angular momentum transport in the solar wind, is explored. Using a 3-fluid model, it is found that 
the radial component of the momentum equation for ions is modified by the force introduced by the 
azimuthal components. This force plays an important role in the force balance in interplanetary space, 
hence impacts the radial flow speed of the species considered, bringing them closer to each other. For 
the fast solar wind, the model cannot account for the decrease of v^p observed by Helios between 0.3 
and 1 AU. However, it is capable of reproducing the profile of Vap measured by Ulysses beyond 2 AU, 
if the right value for Vap is imposed at that distance. In the slow solar wind, on the other hand, the 
effect of solar rotation seems to be more pronounced if one starts with the value measured by Helios 
at 0.3 AU. In this case, solar rotation introduces a relative change of 10-16% in the radial component 
of the fiow speed of the alpha particles between 1 and 4 AU. The model calculations also show that, 
although alpha particles consume only a small fraction of the energy and linear momentum fiuxes of 
protons, they cannot be neglected when considering the proton angular momentum fiux Cp. In most 
examples, it is found that Cp is determined by Vap for both the fast and the slow wind. In the slow 
solar wind, it is also found that the proton and alpha angular momentum fluxes Cp and Ca can be 
several times larger in magnitude than the flux carried by the magnetic stresses Cm- While the sum of 
the angular momentum fluxes Cp = Cp + Ca of both species is found to be smaller than the magnetic 
stress Cm, for the fast and slow wind alike, this result is at variance with the Helios measurements. 
Subject headings: solar wind — Sun: magnetic flelds 



1. INTRODUCTION 

The relative velocity between protons and alpha par- 
ticles in the solar wind, Vap = — Vp, offers important 
clues for understanding the mechanisms responsible for 
the solar wind acceleration. In the fast solar wind with 
proton speeds Vp = \vp\ > 600 km s~^, Helios measure- 
ments made in near-ecliptic regions indicate that there 
exists a substantial Vap = \'Vap\sign(va — Vp) which may 
be 150 km s^^ at the heliocentric distance r « 0.3 AU, 
amounting to ^ 1/4 of the local proton speed. This 



^ap 



decreases with increasing r to 30 — 40 km s 



at 



1 AU. Furthermore, for a continuous high-latitude fast 
stream sampled by Ulysses between May 1995 and Aug 
1996, it was found that the average v^p decreases from 
~ 40 km s-i at 1.5 AU to - 15 km s^^ at 4.2 AU 
(Reisenfeld et al. 2001). As for the slow stream with 



< 400 km s~ 



tends to be zero on average (see 



Marsch et al. (1982) and references therein). However, as 
clearly shown by Figure 11 in Marsch et al. (1982), there 
seem to be two categories of slow solar winds, in one al- 
pha particles tend to flow faster than protons whereas in 
the other this tendency is reversed. That Vap ~ there- 
fore reflects the fact that the two kinds of slow winds 



Electronic address: bblOaber. ac.uk 



have nearly equal opportunities to be present. Marsch 
et al. (1982) interpreted Vap > as a consequence of 
the wave acceleration that favors alpha particles, and 
Vap < as a signature of the solar wind being driven 
primarily by the electrostatic field since alpha particles 
experience only half the electric force that protons do. 
The former interpretation has been corroborated by an 
example on day 117 of 1978, during which period there 
existed simultaneously strong Alfven wave activities as 
well as a significant positive Vap at r w 0.29 AU. This 
Vap (~ 100 km s^"'^) is 20 — 30% of the measured proton 
speed (Marsch et al. 1981). 

Helios measurements have also yielded information 
concerning the angular momentum transport in the so- 
lar wind, in particular the distribution of the angular 
momentum loss C between particles Cp and magnetic 
stresses Cm, and the further partition oi Cp between 
protons Cp and alpha particles Ca (Pizzo et al. 1983; 
Marsch & Richter 1984). In spite of the significant scat- 
ter, the data nevertheless exhibit a distinct trend for 
Cp to be positive (negative) for solar winds with pro- 
ton speeds Vp below (above) 400 km s^^. The precise 
measurement of Ca is even more difficult, however on 
average Ca displays a tendency similar to Cp. The mag- 
netic contribution C]\{, on the other hand, is remarkably 
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constant. A mean value Cm = 1-6 x 10^^ dyne cm sr^^ 
can be quoted for the solar winds of all flow speeds and 
throughout the region from 0.3 to 1 AU. For compari- 
son, the mean values of angular momentum fluxes car- 
ried by ion flows in the slow solar wind are Cp = 19.6 
and £a — 1-3 xlO^^ dyne cm sr^^ (see table II of Pizzo 
et al. (1983)). The fluxes carried by all particles is then 
Cp = Cp + Ca = 20.9 X 10^^ dyne cm sr~^, which tends 
to be larger than Cm- It should be noted that, such a 
value for Cp corresponds to an azimuthal speed of more 
than 7 km for the bulk slow wind, consistent with 
the measurements made before the Helios era (see the 
data compiled in section I of Pizzo et al. (1983)). 

The aforementioned problems are not isolated from 
each other. As a matter of fact, it has been shown that 
in the presence of solar rotation, an additional force ap- 
pears in the meridional momentmn equation of minor ion 
species. The efi'ect of this force is to bring the meridional 
speed of minor ions to that of protons (McKenzie et al. 
1979; Hollweg & Isenberg 1981). It is noteworthy that 
when using the azimuthal speed of minor ions, McKen- 
zie et al. (1979) were mainly concerned with deriving the 
so-called "rotational force" rather than the angular mo- 
mentum transport. Hence the azimuthal proton speed 
was neglected altogether. Hollweg & Isenberg (1981) of- 
fered a more convenient and more self-consistent deriva- 
tion of that force by working in the corotating frame, 
however the azimuthal dynamics was again neglected. 
An extension to the study of McKenzie et al. (1979) has 
been recently given by Li & Li (2006) (hereafter paper 
I) who worked in the inertial frame and, in improvement 
of McKenzie et al. (1979), treated the protons and al- 
pha particles on an equal footing. The resultant model 
is in effect a three-fluid version of the model of Weber 
& Davis (1967). A low-latitude fast solar wind solution 
was worked out, and it was shown that the solar rota- 
tion introduces a barely perceptible difference in the pro- 
files of the meridional ion speeds within 1 AU. However, 
the proton-alpha differential streaming plays a decisive 
role in determining the azimuthal speeds, and thus the 
individual angular momentum fluxes, of not only alpha 
particles but also protons. 

Li & Li's (2006) model was restricted to the region 
within 1 AU, and gave only one solution for the purpose 
of presenting a general analysis on the angular momen- 
tum transport in a three-fluid solar wind. In this study 
we explore, in a quantitative and systematic manner, the 
interplay between the angular momentum transport and 
proton-alpha differential streaming in both the slow and 
fast solar wind for the region extending from the coronal 
base out to 4.5 AU. It should be noted that so far the only 
available measurements of the alpha angular momentum 
flux are from Helios which explored the near ecliptic solar 
wind between 0.3 and 1 AU (Pizzo et al. 1983; Marsch & 
Richter 1984). The Helios measurements will therefore 
be compared with the computed ion angular momentum 
fluxes for the region within 1 AU only. The numerical 
results for the region outside 1 AU will concern only the 
proton-alpha differential streaming Vap, therefore allow- 
ing a comparison of Vap with the Ulysses measurements. 

Before proceeding, we note that several mechanisms 
have been proposed to account for the observed evolu- 
tion of the proton-alpha differential speed Vap in the fast 



solar wind beyond 0.3 AU. Since they have been critically 
reviewed in the introduction section of Kaghashvili et al. 
(2003), some brief remarks are sufficient here. First, the 
pondermotive force due to Alfven waves tends to limit 
\vap\, but it operates too slowly to reduce Vap in the 
observed manner, especially in the region between 0.3 
and 1 AU. Second, various microinstabilities, in partic- 
ular the magnetosonic one, are expected to be opera- 
tional when they have a threshold of Vap of the order 
of the local Alfven speed. However, these microinsta- 
bilities are more likely to be effective in regions beyond 
1 AU where the parallel proton beta is large. For a more 
comprehensive review of these kinetic aspects, please see 
Marsch (2006). Third, the compressional waves accom- 
panying large amplitude Alfven waves may help convert 
the free energy that derives from the differential stream- 
ing into alpha heating, thereby reducing v^p (Kaghashvili 
et al. 2003). Finally, although intuitively appealing, the 
Coulomb friction proves inefficient for regulating Vap in 
the majority of the solar wind measured by Helios. It 
may play some role only for the relatively dense slow 
wind when \vap\ < 15 km s~"'^(sce Figs. 13 and 15 in 
Marsch et al. 1982). The functional dependence on \vap\ 
of the Coulomb collision frequency between protons and 
alpha particles, contained in the coefficient Cq in Equa- 
tion (1), dictates that the friction force will rapidly de- 
crease rather than increase with increasing \vap\ if \vap\ 
exceeds a critical value. Except for the Coulomb friction, 
the aforementioned mechanisms are not included in the 
present study as we try to isolate the effect of solar rota- 
tion on the evolution of Vap, since it is an inherent process 
whose existence is independent of the wave processes. We 
further note that the effect of solar rotation has been in- 
corporated in a number of three-fluid solar wind models 
(e.g., Isenberg & Hollweg (1983), Biirgi & Geiss (1986), 
Hu & Habbal (1999)). However, in all these papers, the 
spiral magnetic field is assumed rather than computed, 
and the azimuthal speeds of ions are neglected. These 
models are therefore not suitable for our purpose, i.e., 
to examine the angular momentum transport in a three- 
fluid solar wind. 

The paper is organized as follows. We flrst give a brief 
overview of the model in section 2. Some further details 
on the numerical implementation of the model are given 
in Section 3. Then section 4 presents the numerical re- 
sults, which are summarized in section 5. 

2. PHYSICAL MODEL 

The solar wind model consists of three species, elec- 
trons (e), protons (p), and alpha particles (a). Each 
species s (s = e,p,a) is characterized by its density 
Tig, velocity v^, mass nis, electric charge e^, and tem- 
perature Ts- The electric charge is also measured in 
units of electron charge e, i.e., = ZsC with = — 1 
by deflnition. The ion mass number Ah follows from 
the relation rrik = AkTUp {k = p,a). The mass den- 
sity of species s is ps = n^mg, and the species par- 
tial pressure is ps — UsksTs, where ks is the Boltz- 
mann constant. From quasi-neutrality it follows that 
He = Hp + ZaUa- We also assume quasi zero-current, 
i.e., Ve = {UpVp + 

Z^Tia^a) /'Ti'ej except when the ion 
momentum equations are derived. The governing equa- 
tions, which self-consistently take into account the effect 
of solar rotation, are identical to those in paper I and 
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will not be given here. They arc appropriate for a time- 
independent solar wind assuming azimuthal symmetry, 
i.e., 9/5^ = in a heliocentric spherical coordinate sys- 
tem (r, 6, <p). In what follows, we will describe briefly 
the model, with the emphasis on the possible effects in- 
troduced by the azimuthal components. In addition, we 
will also give an analysis of the distribution of the total 
angular momentum flux between the magnetic field and 
particles. 

2.1. Model Description 

It was shown in the appendix of paper I that no mass or 
energy exchange between different flux tubes is possible, 
thereby allowing the system of equations to be expressed 
as a force balance condition across the poloidal magnetic 
field coupled with the conservation equations along it. 
In this study, the force balance condition is replaced by 
prescribing a radial (i.e., perfectly monopolar) magnetic 
field. The model equations can therefore be solved for the 
radial distribution of the densities n/j and radial speeds 
Vkr of ion species {k = p,a), the temperatures Tg of all 
species (s = e,p, a), as well as the azimuthal components 
of the magnetic field and ion velocities Vk^,- As such, 
the model can be seen as a 1.5-dimensional (1.5-D) one in 
that the only independent variable is the heliocentric dis- 
tance r, however both the radial and azimuthal compo- 
nents of vectors are retained. To illustrate the effects of 
solar rotation, we will also compute models without the 
(f> components. For the case of description these mod- 
els will be called 1-D ones. Furthermore, the ion heat 
fluxes are neglected, and the field-aligned electron heat 
flux Qe is assumed to follow the Spitzer law, which re- 
sults in V • qe = -{l/a)id/dr)[aKT!^^ (dTe/dr) cos^ $]. 
Here the cross-sectional area of the flux tube a scales as 
a oc 1/Br oc r^, and the magnetic azimuthal angle $ 
is defined by tan$ = B^/Br- As for the electron con- 
ductivity K, the Spitzer value is used, «; = 7.8 x 10"'^ 
erg pjji-i g-i (Spitzer 1962). 

Prom paper I, the equation governing the radial speed 



of ion species 

dVkr 
Vkr^— 

or 



k is 



1 dpk 



Zk dpe GMq 



rikmk dr 
+ flfc + -; — —CQ{vj 



AkUe 



rienik dr r 
Vkr) sec^ $ 



+ 



tan ^Vkr 



dvk(i 
dr 



+ 



r ) 



(1) 



where the subscript j denotes the species other than fc, 
i.e., j = a iov k = p and vice versa, and cq is a co- 
efficient associated with Coulomb frictions. Moreover, 
G is the gravitational constant, Mq is the mass of the 
Sun, and Ofe stands for the acceleration from some ex- 
ternal process. Apart from ak, the meridional accelera- 
tion of the ion flow includes the ion and electron pres- 
sure gradient forces (the terms proportional to dpk/dr 
and dpe/dr^ respectively), the gravitational force, the 
Coulomb friction corrected for the spiral magnetic fleld, 
and the force associated with azimuthal flow speeds (the 
term in the square parentheses, hereafter referred to as 
the "azimuthal force" for brevity) . Note that the electron 
pressure gradient force is part of the electrostatic force 
which is apportioned between two ion species according 
to their charge-to-mass ratios Zk/Ak {k =p,a). 



Introducing (j) components impacts the meridional dy- 
namics through two means. First, the electron conduc- 
tivity is in effect reduced by a factor cos^ $, thereby af- 
fecting the electron temperature and hence the electron 
pressure gradient force. Second, the azimuthal force en- 
ters into the meridional momentum equations. The two 
effects can both be important, as will be illustrated by 
the numerical solutions presented in Section 4.2. 

2.2. Angular Momentum Loss Rate and Its Distribution 
between Particles and Magnetic Stresses 

In what follows, a simple analysis based on constants 
of motion for the governing equations shows how the an- 
gular momentum loss rate per steradian C is distributed 
between particles C p and magnetic stresses Cm in a solar 
wind where different ion species flow with different veloc- 
ities. As such, this discussion complements that of Weber 
& Davis (1967) and Marsch & Richter (1984) where the 
solar wind is seen as a bulk flow. 

Following paper I, the expressions for the azimuthal 
components Vpcf,, v^^ and B^ are 



fir sin 



Va<t> = 




Ml 



Ba = Qr sin 9 



4TTppVpr{l + 0r\/r^ -1 



Br 



(2c) 



in which v^p^r is the radial component of the velocity 
difference vector v^p, 17 is the angular rotation rate of the 
flux tube, and the constant ( = (paVar) / (PpVpr) denotes 
the ion mass flux ratio. The terms designated by II are 
associated with the differential streaming, whereas the 
terms denoted by I take care of the rest. The subscript 
A denotes the Alfvenic point where Mt = 1, Mt being 
the combined meridional Alfvenic Mach number deflned 
by 



Ml 



Ml 



B^/Anpk' 



(3) 



with k = p,a. 

The field and particle contributions to the angular mo- 
mentum loss rate per steradian C are 



Cm = —r sin9- 



,B^Br 

47r ' 



Ck = r^ sin 0pkVkrVk4>, (4) 



where k =p,a. Prom Equation (2) follows that 

C = MQ.r\sm.e\, (5) 

in which M = (1 + QppVprr'^ is the total mass loss rate of 
the solar wind. The ratio Ckj Cm is rather complex and 
had better be examined case by case due to the presence 
of terms labeled II in Equations (2a) and (2b). Note that 
terms II derive essentially from the requirement that the 
proton-alpha velocity difference vector be aligned with 
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the instantaneous magnetic field (Li & Li 2006). The az- 
imuthal speeds are expected to be determined by these 
terms when there exists a substantial Vap,r, which is al- 
most exclusively determined by the way in which the ex- 
ternal energy is distributed among different species in the 
corona. As shown in section 4, such a situation happens 
for the majority of the numerical solutions in interplane- 
tary space. Hence the azimuthal speeds and the resultant 
specific angular momenta of ion species measured in situ 
are essentially determined by the processes that happen 
in the corona. Furthermore, it follows from the identity 
Pp'i'irrM^^/'L'ar = Pai'ar / Vp,. that tcrins II havc no con- 
tribution to the overall angular momentum loss rate. As 
a result, a concise result can be found for the ratio of 
the overall particle contribution Cp = Cp + Ca to the 
magnetic one Cm- By defining the bulk velocity 

u=(vp + Cv„)/(l + C), (6) 
one can see M|. = {r/rA)^{ur/ur^A)- As a result, 

Cp _ Ur/Ur,A - 1 ,„s 



1 - r\lr^ 



Vp, the results are identical to those in Marsch & 



It is interesting to see that only one single speed Ur enters 
into the ratio Cp/Cm- 
In comparison with Weber fc Davis (1967) and Marsch 

& Richter (1984), the expressions here indicate that, a 
proper definition of the Alfvenic point is through Mt = 
1, and the bulk flow velocity u should be defined as an av- 
erage between the proton and alpha velocities, weighted 
by their respective mass fiuxes (cf. Eq.(6)). Of course 
these differences only matter when the abundance of al- 
pha particles and the differential streaming between the 
protons and alpha particles are large. When ^ = or 
v< 

Richter (1984) 

3. NUMERICAL IMPLEMENTATION 

The model has been briefly described in the preced- 
ing section where a rather general analysis is also given. 
However, it is still necessary to solve the governing equa- 
tions numerically for a quantitative analysis to be made. 
In particular, some external heating and/or momentum 
addition need to be applied to the ions to generate a 
solar wind solution. The implementation of the energy 
deposition and a description of the method of solution 
are given in this section. 

3.1. Energy Deposition 

The external energy deposition is assumed to come 

from an ad hoc energy flux which dissipates at a constant 
length Ici- The resultant total heating rate is therefore 



Br 



■ exp 



I, 



where Fe is the input flux scaled to the Earth orbit Re = 
1 AU, BrE = 3.37 is the radial magnetic field strength at 

Re, and Rq is the solar radius. The energy is deposited 
in protons entirely as heat Qp, but deposited in the alpha 
gas in the form of both heating Qa and acceleration Oq,. 
To be more specific, Q is apportioned among Qp, Qa and 
tta in the following way. 



Qp ~l~ Qo 



Q, — = 
Qp 



XPa 
Pp 

Qa 



TABLE 1 

Parameters characterizing the energy 

dkrositton to ton ^pecie^ 



wind 


Fe 


h 


X 


Xd 


type 


(erg cm ^ s ^) 


(Re) 






Fast 


1.8 


5 


1.8 


see note'' 


Slow 


1.2 


1 


2.9 






^Xd decreases with r from 30 to 1 with a sharp 
transition at 2 Rq 



Xd 
PaVo' 



(8b) 



where vq = yHtBTp/rnp is some characteristic speed. 

The choice of the heat and momentum deposition needs 
some explanation. The exponential form of Q was first 
suggested by Holzer & Axford (1970) and was later em- 
ployed in a large number of studies. The form of Q 
adopted here is slightly different from the original ver- 
sion to ensure that it mimics the dissipation of a flux of 
non-thermal energy, for instance, the dissipation of low- 
frequency Alfven waves. This may happen as a result 
of a turbulent cascade towards high parallel wavenum- 
bers where the wave energy is picked up by ions through 
the ion-cyclotron resonance. (For more details, please 
see the review by HoUweg & Isenberg (2002) dedicated 
to this topic.) Once Q is specified, the way it is dis- 
tributed among different ion species, characterized by 
the parameter x, depends on the resonant interaction 
between cyclotron waves and ions. Note that x indicates 
how alpha particles are favored when Q is distributed, 
with X = 1 standing for the neutral heating: the total 
energy that goes to ion species k is proportional to its 
mass density pk [k = p,a). Previous computations in- 
volving ion cyclotron waves indicate that, in the case of 
neutral heating, the alpha particles tend to flow slower 
than protons (see the dispersionless case in Fig.l of Hu 
& Habbal 1999). Only when the alphas are energetically 
favored in the corona can the modeled Vap be positive in 
interplanetary space. This happens when x > 1- (Note 
that we adopted a constant x throughout the computa- 
tional domain for simplicity.) Listed in Table 1, and 
unless otherwise stated, the heating parameters are cho- 
sen to produce fast or slow solar wind solutions with 
realistic ion mass fluxes and terminal speeds. 

3.2. Method of Solution and, Boundary Conditions 

The governing equations (Eqs.(l) to (7) in paper I) are 
cast in a time-dependent form and are solved by using a 
fully implicit numerical scheme (Hu et al. 1997). Start- 
ing from an arbitrary guess, the equations are advanced 
in time until a steady state is reached. The computa- 
tional domain extends from the coronal base (1 Rq) to 
4.5 AU. At the base, the ion densities and species tem- 
peratures are fixed, ng = 3 x 10* cm^'', Ua/up = 0.1, 
Te = Tp = Ta = 1.2 X 10^ K. The radial components of 
ion velocities, Vpr and Var, are specified to ensure mass 
conservation. On the other hand, Vp^ and are spec- 
ified according to Equation (2). At the outer boundary 
(4.5 AU), all dependent variables are linearly extrapo- 
lated for simplicity. We take f2 = 2.865 x 10~^ rad s~^, 
which corresponds to a sidereal rotation period of 25.38 
days. 

4. NUMERICAL RESULTS 
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As has been described in the introduction, we are in- 
terested in answering two questions: 1) To what extent 
is Vap, the differential streaming between protons and 
alpha particles, affected by introducing solar rotation? 
and 2) Is it possible to reconcile the model results with 
the Helios measurements, as far as the specific angular 
momentum fluxes are concerned? In this section, we will 
first examine the fast solar wind and then move on to 
the slow one. 

4.1. Fast Solar Wind Solutions 

4.1.1. Effect of Solar Rotation on the Differential 
Streaming 

Given in Figure 1 is a comparison of the 1.5-D (solid 
curves) with the 1-D model (dotted curves) of the fast 

solar wind solution in the equatorial plane (colatitude 
9 = 90°). Plotted are the radial profiles for (a) the 
radial component of the proton-alpha velocity differ- 
ence Vap,n and (b) the radial ion speeds Vpr and Var- 
In Fig. la, the dash-dotted line shows the parameter 
Vap = UQ,p,j.sec$, which also takes into account the az- 
imuthal component of the difference vector. The asterisk 
denotes the Alfvenic point in the 1.5-D model, which lies 
at ta = 12.6i?0. For comparison, the open boxes give 
the near-ecliptic Helios measurements of w„p for the fast 
wind as reported by Marsch et al. (1982). 

From Figs. la and lb, one can see that introducing 
the solar rotation only affects the solar wind beyond 
say 50 Rq. Thus it is not surprising that the ion 
fluxes are unchanged in the 1.5-D model, as compared 
with the 1-D one. Both models yield an ion flux ra- 
tio of {uaVar) / {npVpr) = 0.026, and a proton flux of 
{npVpr)E = 2.72 X 10® cm~^ s~^ when scaled to 1 AU. 

It can be seen from Fig.la that, although Vap,r de- 
creases with radial distance r for both the 1-D and 1.5-D 
models beyond their mutual maximum of 147 km s~^ 
attained at 9.4 Rq, the reduction in Vap,r is more promi- 
nent in the 1.5-D model from 0.3 to 4 AU. The reduction 
in Vap for the 1.5-D model, when compared with the 1-D 
one, is also substantial, although less significant than the 
reduction in Vap.r due to the presence of the azimuthal 
component of v^p. Nevertheless the 1.5-D model yields 
a Vap of 11.4 km s~^ at 4 AU, as opposed to 29.3 km s~^ 
attained in the 1-D model. From Fig. lb one can see that 
the reduction in Vap,r in the 1.5-D model is achieved by 
lowering the radial speed profile of alpha particles Var 
and raising that of protons Vpr- This effect is barely 
perceptible at 0.3 AU but becomes more obvious with 
increasing r. 

Figure la indicates that the modeled Vap profile de- 
viates significantly from the Helios measurements, espe- 
cially in the region closer to the Sun. In other words, 
the deceleration of alpha particles relative to protons as 
measured by Helios cannot be explained by merely invok- 
ing the (j) components. However can this mechanism be 
entirely ruled out as far as the Ulysses measurements are 
concerned? Such a possibility is examined in Figure 2. 
In Fig. 2a, the distribution of Vap in the r — 6 plane for 
the region extending from 1 to 4.2 AU is displayed as 
contours equally spaced by 1 km s~^. The solutions are 
obtained by varying the colatitude alone. Recalling 
that by assumption the solar wind flows in a perfectly 
monopolar magnetic field, one may see that the latitu- 



dinal dependence is entirely due to solar rotation. From 
Fig. 2a it is obvious that in the regions close to the pole, 
say < 10°, Vap shows little radial dependence. How- 
ever with increasing 6, the radial gradient in Vap becomes 
increasingly significant. 

In Fig. 2b, the solid curve represents the modeled Vap 
profile along the Ulysses trajectory delineated by the 
dashed curve in Fig. 2a in the interval between 1995 May 
05 and 1996 Aug 02. During this period, Ulysses sam- 
pled a continuous, undisturbed, high speed stream above 
the latitude of 30°N in the radial range of 1.5 ~ 4.2 AU. 
The actual Ulysses measurements are given by the error 
bars which correspond to the 25-75% percentiles of the 
Vap distribution over 10-day bins (Reisenfeld et al. 2001). 
It can be seen that, for r < 2 AU, Vap in the numerical 
solutions increases as opposed to the observed tendency 
for Vap to decrease. This is because during this period, 
Ulysses was moving towards the pole and the effect of so- 
lar rotation diminishes with decreasing 9. However, for 
r > 2 AU, the decrease in Vap from 30.7 km s~^ at 2 AU 
to 12.4 km s~^ at 4.2 AU is consistent with the observed 
values within the accuracy of the measurements. From 
this we conclude that solar rotation should not be ne- 
glected in attempts to understand the Ulysses measure- 
ments of Vap beyond 2 AU. It is necessary to stress that 
this conclusion should not be confused with that reached 
by Reisenfeld et al. (2001) who demonstrated by showing 
a profile designated "Ulysses rotational deceleration" in 
their Figure 8 that the Vap profile measured by Ulysses 
has little to do with the rotational deceleration. In fact, 
we agree with these authors that introducing solar rota- 
tion cannot provide a unified mechanism to account for 
the observed deceleration of Vap from 0.3 AU onwards. 
What Fig. 2b shows is that, if a reasonable Vap is given at 
2 AU, the solar rotation alone can decelerate the alpha 
particles relative to protons in the observed fashion. 

4.1.2. Angular Morn,6nturn, Transport 

We now turn to the problem of the angular momentum 
transport in the fast solar wind along 9 = 90°. Figure 3 
describes the radial distribution of (a) the azimuthal ion 
speeds Vp^f, and Va<p, and (b) the specific angular mo- 
mentum fluxes including those carried by protons Cp, 
alpha particles Ca and magnetic stresses Cm- The to- 
tal angular momentum flux C should be a constant, as 
is indeed the case in the numerical solution. To isolate 
the contribution of the differential streaming to individ- 
ual azimuthal speeds, term I (II) in Eqs.(2a) and (2b) 
is plotted as the dashed (dash-dotted) curve between 0.3 
and 4 AU. In Fig. 3b, the dashed line is used for plotting 
negative values. The asterisks in both panels mark the 
Alfvenic point va- 

From Fig. 3a one can see that with the development of 
Vap,r, some substantial difference arises between Vp^ and 
Va(j>- To be more specific, Vp^f, increases from the coro- 
nal base to 5.84 km s~^ at 6.9 Rq, beyond which Vp^ 
decreases gradually to 1.1 km around 4 AU. On the 
other hand, Vacf, reaches a local maximum of 4.3 km s^^ 
at 2.88 Rq. Rather than partially corotate with the Sim, 
alpha particles then gradually develop an azimuthal ve- 
locity in the direction of counter-rotation with the Sun: 
is negative beyond 11.7 Rq. Then Va,p attains a lo- 
cal minimum of -14.4 km s~^ at 324 Rq and increases 
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Fig. 1. — Effect of solar rotation on a 3-fluid fast solar wind in the equatorial plane d = 90°. The aziinuthal components of ion velocities 
and the magnetic field are incorporated self-consistently in the 1.5-D model (solid curves), but excluded in the 1-D model (dotted curves). 
Radial profiles are shown for several flow parameters, (a) Radial component of the proton-alpha relative velocity Vcp = Vc — Vp, and (b) 
the radial speeds of protons Vpr and alpha particles Var- In panel (a), Vap = Va-p,r sec# is given by the dash-dotted line, where <I> is the 
magnetic azimuthal angle, deiined by tan<I> = B^/Br- The asterisk refers to the Alfvenic point where Mj- = 1, My being the combined 
Alfvenic Mach number defined by Equation (3). In addition, the open boxes give the near-ecliptic Helios measurements of Vap reported by 
Marsch et al. (1982). 
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Fig. 2. — Effect of solar rotation on a 3-fluid fast solar wind, (a) Distribution in the r — 6 plane of the proton-alpha differential speed 
Vap for the region extending from 1 to 4.2 AU. The contours are equally spaced by 1 km s~^. The dashed line depicts the trajectory of 
Ulysses in the interval between 1995 May 05 and 1996 Aug 02. The solutions are obtained by varying the colatitude 6 alone with all other 
parameters unchanged. As such, the latitudinal dependence of Vap is entirely due to solar rotation, (b) Radial proflle of the modeled 
Vap along the Ulysses trajectory (solid line). The error bars represent the 25-75% percentiles of the Vap distribution over 10-day bins as 
sampled by Ulysses (Rciscnfcld ct al. 2001). 



thereafter to -9.35 km at 4 AU. It can be seen that, 
for r > 0.3 AU not only Va(t> but also Vp^ are determined 
mainly by term II (the dash-dotted lines) in their respec- 
tive expressions (2a) and (2b). 

Despite the small values of Up^, Fig. 3b reveals that 
Cp exceeds Cm in the region beyond 129 Rq. Further- 
more, Ca becomes larger than Cm in magnitude for 
r > 186 Rq. Consider the values at 1 AU. As part 
of the total flux C = 2.47 (all values are in units of 
10^'' dyne cm sr^^), the magnetic part Cm = 1-72 is 
smaller than the proton contribution Cp — 2.86 which 
however is largely ofi'set by the alpha contribution Ca = 
—2.12. As a matter of fact, magnetic stresses are al- 
ways the primary constituent in C, the particle con- 
tribution clS £i whole is smaller than Cm throughout 
the computational domain. At 1 Rq, 99.37% of C is 
contained in magnetic stresses Cm- Although Cm de- 
creases monotonically with increasing distance, it still 
amounts to 69.2% of £ at 4 AU. This is understand- 
able in view of Equation (7). For r ^ Rq, one can 
see that r <g; and Ur <^ Ur a- It then follows that 



CM/C^l-{r/rA) 
one can see Cm/C 



i 1. On the other hand, for r ^ ta 
Ur^A/ur- That Cm takes up most 
of the total flux C asymptotically therefore stems from 
the fact that Ur^A/ur varies little beyond ta- 

The value oi Cm (in units of 10^^ dyne cm sr^^) shows 
only a slight decrease from 1.76 at 0.3 AU to 1.72 at 1 AU. 
This is consistent with the quoted value of 1.6 ~ 1.9 as 
measured by Helios (Pizzo ct al. 1983; Marsch & Richter 
1984). However, Helios measurements indicate that in 
the fast solar wind, both Ca and Cp tend to be negative. 
As a result, the overall particle contribution £p is neg- 
ative, which can be realized only if the solar wind bulk 
speed Ur decreases with r in the region r > va- This is 
not possible by simply varying the hc;ating parameters 
in the model. Instead, we would rather intcrpr(^t Fig. 3 
as an indication that the the proton angular momentum 
flux cannot be studied on its own, i.e., by simply incor- 
porating (f) components into models of an electron-proton 
plasma. Although alpha particles take up only a minor 
part of the solar wind momentum and energy fluxes, the 
proton-alpha diS'erential speed is the primary contributor 
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Fig. 3. — Angular momentum transport in a 3-fluid fast solar 
wind, (a) Azimuthal speeds of protons v^^ and alpha particles 1)^0 1 
and (b) specific angular momentum fluxes carried by protons yCp, 
alpha particles Ca and magnetic stresses Cm ■ The fluxes add up to 
C In panel (a), the dashed line shows the contribution from term 
I in Eqs.(2a) and (2b), whereas the dash-dotted curves give the 
contribution from term II. In panel (b), the dashed line represents 
negative values. The asterisks mark the position of the Alfvenic 
point TA- 

to individual azimuthal ion speeds and therefore individ- 
ual ion angular momentum fluxes. 

4.2. Slow Solar Wind Solutions 

The effect of solar rotation on individual ion speeds in 
the fast solar wind is rather weak inside 1 AU. Is this 
also the case for the slow wind? Moreover, the Helios 
measurements show that in the slow solar wind, the pro- 
ton and alpha angular momentum fluxes are both posi- 
tive. When added together, they are also larger than that 
carried by magnetic stresses. Can this measurement be 
explained by the present model? To answer these ques- 
tions, we will first examine an example with a significant 
positive Vap^r to gain some insight. A parameter study 
surveying the parameter x will then be presented. 

4.2.1. Effect of Solar Rotation on the Differential 

Streaming 

Figure 4 compares a 1.5-D (solid curves) with 1-D 
model (dotted curves) of the slow solar wind in the equa- 
torial plane 9 = 90°. Plotted are the radial profiles for 
(a) the radial component of the proton-alpha velocity 
difference v^p.r, and (b) the radial speeds of ions v^.^ 
{k = p,a). In Fig. 4a, the asterisk denotes the Alfvenic 
point in the 1.5-D model, which now lies at rA = 14.2 Rq. 
The dash-dotted line is used to show the parameter Vap- 

As in the case for the fast wind, introducing the so- 
lar rotation does not alter the ion fluxes. When scaled 
to lAU, both the 1-D and 1.5-D models yield a proton 
flux of {npVpr)E = 3.36 x 10^ cm^^ s^^, and a flux ratio 
of {uaVar) I {npVpr) = 0.024. Furthermore, introducing 
the (!> components reduces the proton-alpha differential 
speed by decelerating the alphas and accelerating the 
protons. Take the radial speeds at 1 AU for instance. 
The 1.5-D model yields Var = 389 km s~^, whereas the 



1-D model results in Var — 435 km s ^ 
1.5-D (1-D) model yields 353 (347) km s'^ As a result, 
Vap,r at 1 AU is 88 km s~^ for the 1-D but 36 km s~^ 
for the 1.5-D model. Going further to 4 AU, in the 1-D 
model alpha particles still flow considerably faster than 
protons: values of v^r = 435 and Vpr = 353 km s^^ are 
found at 4 AU. In contrast, Var and Vpr attain nearly 
the same value (~ 365 km s~^) in the 1.5-D model. It 
should be noted that, a value of Vap = 88 km s~^ at 
0.3 AU seen in Fig. 4a for the 1.5-D model is not un- 
realistic. Even larger values have been found for Vap 
by Helios 2 when approaching perihelion (Marsch et al. 
1981). However Vap at 1 AU seems somehow tuirealis- 
tic: A value of 56 km s~^ is substantially larger than the 
local Alfven speed Va = Br seci>/ ^/4Tr{pp + pa) which 
yields 35 km s~^. By contrast, the Helios measurements 
indicate that in the slow wind the ratio v„p/V^ spans 
a broad distribution between —1 and 1 but peaks at 
~ ±0.25 (Figure 11 in Marsch et al. (1982)). The discrep- 
ancy indicates once again that invoking the solar rotation 
alone is not an efficient way to decelerate the alpha par- 
ticles relative to protons. However, if some yet unknown 
mechanism results in a substantial Vap at 0.3 AU as ob- 
served, the effect of solar rotation on the evolution of 
Vap with radial distance r beyond 0.3 AU probably can- 
not be neglected from the perspective of slow solar wind 
modeling. For instance, in this particular solution, the 
changes introduced in the alpha speed Va,- at the distance 
of 4 (1) AU can be 16.1% (10.5%) the values in the case 
without taking solar rotation into account. 

Returning to the discussion of Vcep, one may expect 
that Vap is largely determined by the proportion of the 
external energy deposited to alpha particles, character- 
ized by the parameter x- Figure 5 displays the depen- 
dence on X of Vap at several different radial distances 
as labeled. Results from 1.5-D models are displayed by 
thick curves, whereas those from 1-D models are given by 
thin curves for comparison. Inspection of the solid lines 
shows that the transition from negative to positive Vap 
at 0.3 AU occurs at x = 1.6 in both models. With x var- 
ied between 0.5 and 3.4, Vap at 0.3 AU goes from —132 
to 127 km s"^ for the 1.5-D model, while from -143 to 
134 km s"'^ for the 1-D model. When comparing the thin 
and thick solid lines, one can see that the profile of Vap 
for 1-D models differs only slightly from that for 1.5-D 
ones. Going from 0.3 to 4 AU, the magnitude of Vap is 
significantly reduced by introducing the azimuthal com- 
ponents. Unsurprisingly this effect is more prominent 
at the high and low ends. For instance the 1-D model 



for X = 3.4 yields Vo 



120 km 



at 4 AU, whereas 



the corresponding 1.5-D model gives Vap = 21.8 km s~^. 
It turns out that introducing the components reduces 
Vap via different mechanisms for different portions of this 
parameter range. For x between 1.0 and 2.3, the proton- 
alpha friction eventually becomes operative between 0.3 
and 4 AU to suppress any significant \vap\- As a re- 
sult, no obvious deviation shows up between 1-D and 
1.5-D models. If x > 2.3, the reduction of \vap\ in the 
1.5-D model is achieved partly by the enhanced electron 
pressure gradient force (which stems from the enhanced 
electron temperature due to reduced effective electron 
thermal conductivity in presence of the spiral field) and 
more importantly via the azimuthal force. On the other 
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Fig. 4. — Effect of solar rotation on 3-fluid slow solar winds in the equatorial plane 9 = 90°. The 1.5-D model is given by solid curves 
whereas the 1-D one is given by dotted curves, (a) Radial component of the proton-alpha relative velocity v^p and (b) The radial speeds of 



protons Vpr and alpha particles Var- In panel (a), the asterisk refers to the Alfvenic point, and the dash-dotted line shows Vap ■ 
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Fig. 5. — Dependence of the proton-alpha differential speed Vap 
on the parameter x at three radial distances, as labeled, for the 
1.5-D model (thick curves) and the 1-D one (thin curves). 
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Flc;. 6. — Values of several parameters at 1 AU as a function 
of parameter x- (a) Azimutlial speeds of protons Vp^ and alpha 
particles v^^, and (b) Specific angular momentum fluxes carried 
by the proton fluid Cp, the alpha fluid and magnetic stresses 
Cm- The total flux C is also given. 



hand, for x < 1, the reduction in \vap\ is solely due to the 
azimuthal force, since in this case, the electron pressure 
gradient force tends to increase the magnitude of "Vfj^p clS ct 
consequence of the smaller charge-to-mass ratio of alpha 
particles compared with that of protons. 

4.2.2. Angular Momentum Transport 



Figure 6 displays the distribution with x of the val- 
ues at 1 AU of (a) the azimuthal speeds of ions Vpcj, and 
Vad), and (b) specific angular momentum fluxes includ- 

alpha particles Ca, as 
M- From Fig. 6a one can 
at 1 AU tend to have opposite 



ing that carried by protons Cp 
well as by magnetic stresses C 



see that 



and 



signs, an evidence that they are mainly determined by 
the terms associated with Vap,r in Eqs.(2a) and (2b). 
As a matter of fact, for all x but x = l-Q? term II in 
Equation (2b) contributes more than term I to Vact>- As 
to Vp,p, with the exception of 1.4 < x < 2, the main 
contribution also comes from term II in Equation (2a). 
Looking at Fig. 5, one finds that for this narrow range of 
X, Vap at 1 AU ranges from -3 to 1.1 km s~^. At substan- 
tial relative speeds between protons and alpha particles, 
Vp^ can become as large as 5.84 km s~^ when x = 3.4, 
for which Va^ = —58.5 km s~^. At the other extreme 
X = 0.5, Vpc/, (va^) is -4.1 (59.3) km s~^. In addition, 
from Fig. 6b one can see that Cm or C varies little, al- 
though Cp and Ca vary significantly within the parame- 
ter range as a consequence of the sensitive dependence of 
Vatj, and Vp^ on x- For instance, when x varies from 0.5 
to 3.4, C shows only a slight increase from 3.39 to 3.86 
xlO^^ dyne cm sr"-*^, and Cm increases from 3.09 to 3.51 
xlO^^ dyne cm sr^^. The increase of C with x can be 
explained by Equation (5). For the solutions considered, 
the location of the Alfvenic point varies little. However, 
the mass loss rate, which is mainly determined by the 
proton flux, increases with increasing x as a consequence 
of more energy being deposited to the subsonic portion 
for the proton gas (Leer & Holzer 1980). On the other 
hand, the ratio Cm/C shows little variation with x- This 
can be seen as an indication that the ratio of the bulk 
speed Ur attained asymptotically to that at the Alfvenic 
point is rather insensitive to x, since Cm/C ~ Ur^A/ur 
in the region r ^ ta (cf. Eq.(7)). Moreover, if m,. varies 
little beyond r^, as is the case for all the obtained slow 
wind solutions. Cm should be only slightly larger than 
unity for r ^ ta- The expectation is reproduced by 
Fig. 6b from which one can see that magnetic stresses are 
always the most important term in C, despite the fact 
that individual angular momentum fluxes carried by ion 
species, Cp or Ca, can be much larger than Cm- 

Now the model results can be compared with the He- 
lios measurements for specific angular momentum fluxes 
(Pizzo et al. 1983; Marsch & Richter 1984). One can see 
that, although Cm falls within the uncertainties of the 
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measurements, and values like 1.3x10 dyne cm sr~ 
can be found for Ca at x « 1, no values as large as 
19.6x10^^ dyne cm sr~^ can be found for Cp unless the 
alpha particles are heated even more intensely. Moreover, 
Fig. 6b does not reproduce the observed tendency for the 
particle contribution £p to be larger than the magnetic 
one Cm at 1 AU. As has been discussed, the reason is that 
the solar wind bulk flow Ur experiences little acceleration 
outside the Alfvcnic point ■ It is unlikely that this fea- 
ture can be changed by simply varying other model pa- 
rameters. Actually the discrepancy between model pre- 
dictions and measurements on the relative importance of 
Cp and Cm is common to all existing models that treat 
the meridional and azimuthal dynamics sclf-consistently 
(cf. sections I and VII of Pizzo et al. (1983)). To ex- 
amine the possible sources leading to this discrepancy, in 
what follows we consider some of the physics that is not 
accounted for in the present model. 

In the framework of steady state models, the fast 
and slow magnetosonic waves, in the WKB approxima- 
tion, have non-zero rcj) components in their stress ten- 
sors, thereby also contributing to the angular momen- 
tum loss of the solar wind (Marsch 1986). On the one 
hand, this contribution may be comparable to that due 
to the background flow and magnetic field stress when 
the waves have sufficiently large amplitudes. On the 
other hand, depending on the propagation angle with 
respect to the background magnetic field, the particle 
part may dominate the magnetic one in the wave con- 
tribution. (This can be seen if one evaluates the rcj) 
component of Eq.(7) by using, e.g., the fast wave eigen- 
function Eq.(lO) in Marsch (1986)). Therefore, including 
compressional waves may provide a possible means to re- 
solve the apparent discrepancy. Furthermore, extending 
the WKB analysis to higher order, HoUweg (1973) found 
that Alfven waves with finite wavelength may carry a 
non-vanishing flux of angular momentum, but the con- 
tribution seems to worsen rather than improve the dis- 
crepancy. Finally, we refer the reader to the discussion 
of Hu et al. (2003), whose 2.5-D model also incorporates 
the momentum deposition by Alfven waves. Instead of 
attributing this discrepancy to mechanisms missing in 
the steady state model, Hu et al. (2003) proposed that 
the tendency for £p to be larger than Cm is due to inter- 
planetary dynamical processes such as the forward fast 
shocks traversing the solar wind. Their discussion also 
applies to the present model, although the details of the 
response of the three-fluid plasma to fast shocks need to 
be explored. 

5. SUMMARY 

The near-ecliptic measurements of the fast solar 
wind by Helios demonstrated the existence of a signif- 
icant proton- alpha differential speed, Uap, varying from 
150 km s-i at 0.3 AU to 30-40 km s'^ at 1 AU(Marsch 
et al. 1982). A steady decrease of Vap between 1.5 and 
4 AU was also measured by Ulysses out of the ecliptic 
plane (Reiscnfeld et al. 2001). The Helios measurements 
also indicated the existence of two categories of slow so- 
lar wind. In one category, the alpha particles tended to 
be slower than the protons, whereas the opposite trend 
was found in the other category (Marsch et al. 1982). The 
measurements also yielded information about the angu- 
lar momentum transport, in particular, the specific an- 



gular momentum fluxes in the magnetic stress Cm, and 
in the two major ion species Cp and £„; (Pizzo et al. 1983; 
Marsch & Richter 1984). 

In an attempt to account for these observations, we in- 
troduced solar rotation self-consistently in a three-fluid 
solar wind model. The impetus for this approach was 
inspired by the fact that the solar rotation is an inherent 
process, independent of different mechanisms invoked for 
accelerating the solar wind. By exploring the resulting 
interplay between v^p and the angular momentum trans- 
port we reached the following conclusions: 

1. The force introduced in the radial momentum 
equations by the azimuthal components can play 
a significant role in the force balance in interplan- 
etary space. Its main effect is to reduce the differ- 
ence between the ion ffow speeds. 

2. While the effect of solar rotation could not account 
for the decrease in Vap in the fast solar wind from 
0.3 to 1 AU, it could account for the Vap profile 

measured by Ulysses beyond 2 AU if an appropriate 
value for Vap was chosen at 2 AU. 

3. The effect of solar rotation was found to be more 
pronounced for the slow solar wind if a significant 

Vap developed by 0.3 AU, resulting in a relative 
change of 10-16% in the radial speed of the alpha 
particles between 1 and 4 AU. 

4. For the slow solar wind solutions considered, the 

angular momentum flux carried by alpha particles, 
£q,, at 1 AU, was almost exclusively determined by 
Vap- Moreover, Cp was also determined by Vap for 



5. In the fast solar wind, Cp and Cq beyond 0.3 AU 
are also mainly due to Vap- This suggests that even 
though alpha particles take up only a small fraction 
of the energy and linear momentum fluxes of pro- 
tons, they cannot be neglected when the proton an- 
gular momentum flux is concerned, i.e. Cp cannot 
be studied by simply incorporating the azimuthal 
components into models of an electron-proton solar 
wind. 

6. For the slow solar wind, Cp and Ca can be several 
times larger in magnitude than the flux carried by 
the magnetic stresses Cm- However, the total par- 
ticle distribution Cp = Cp + Ca is always smaller 
than Cm for all solutions. While this tendency is 
inherent to the present model, it is at variance with 
the Helios measurements. 

In closing, some words of caution seem in order. As 
has been discussed, a pure fluid model of a quiet solar 

wind as adopted in this study is rather idealized. In 
reality, the solar wind will be subject to various distur- 
bances such as compressional magnetosonic waves, shock 
waves, and other processes resulting from stream-stream 
interactions. These disturbances may significantly mod- 
ify the distribution of the angular momentum flux be- 
tween particles and magnetic field stresses, and may al- 
ter the overall angular momentum loss rate of the solar 
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wind. With increasing distance, such as in the outer he- 
liosphere, these disturbances are expected to become in- 
creasingly important. With this caveat in mind, the pre- 
sented study nevertheless provides a better understand- 
ing of the underlying physics. Even when the above- 
mentioned dynamical processes are taken into account, 
one would expect to see that solar rotation has significant 
effects on the alpha flow speeds in interplanetary space. 



and in turn, alpha particles are important in the problem 
of angular momentum transport of the solar wind. 
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